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components that fulfils system requirements with the lowest system Total Cost of Ownership (TCO). The 
passive wind system associated with the storage bank interacts with wind speed and load cycles. A set of 
small power passive wind turbines spread on a convenient power range (2-16 kW) are obtained through 
Keywords: an IOD process at the device level detailed in previous papers. The system cost model is based on data 
Passive wind turbine sheets for the wind turbines and related to battery cycles for the storage bank. From the range of wind 
A e turbines, a “system level” optimization problem is stated and solved using an exhaustive search. The 

optimization results are finally exposed and discussed through a sensitivity analysis in order to extract 


Battery cycles ; ; Ea A 
Total cost of ownership the most robust solution versus environmental data variations among a set of good solutions. 


Sensitivity analysis © 2013 Elsevier Ltd. All rights reserved. 
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Nomenclature 

Č scale factor (m/s) 

Co cost of one deep battery cycle (€) 

C3 battery element nominal capacity (Ah) 
Cpar battery bank cost (€) 

Ccell battery cell cost (€) 

Cr cycle to failure 

Cwr wind turbine cost (k€) 

DOD depth of discharge 

I, battery element nominal discharge current (A) 
Teej battery cell current (A) 

Ich_max battery cell maximum charge current (A) 


Taisch_max battery cell maximum discharge current (A) 


k shape factor 
N number of cycles at a given DOD 
Neel_p number of cells in parallel 


1. Introduction 


According to the World Energy Agency [1], some 1.5 billion 
people had no access to electricity in the world by 2009, with 
more than 80% of habitants in rural zones. Providing consumers in 
remote areas with reliable and cheap electricity becomes a priority 
in several developed and undeveloped countries such as in the 
case of isolated cities in Tunisia. The steadily increasing demand of 
fossil fuels along with concerns about global warming, presents 
natural renewable energy sources as attractive solutions. Among 
these sources wind energy systems (free in their availability, 
renewable and non-polluting) with storage are among the most 
competitive alternatives for electrifying remote consumers and 
they are widely used in both autonomous or grid connected 
applications. These systems can also operate in parallel with 
others available energy sources (fuel cells, diesel generators) and 
several means of storage (accumulators, H2 storage, etc.) in order 
to enhance the system reliability [2-5]. 

However, the drawbacks of such sources are their TCO which is 
still expensive, especially for small wind turbine systems. More- 
over to assure the service continuity and to protect the battery 
against deep discharges (subsequently extend the battery bank 
life), such systems require an additional dynamic source of energy 
or an optimal wind system design [8,9]). Recently, several 
researches based on global optimization techniques have been 
focused on the design of optimal system configurations which 
meet the load demand for given climatic data [10-12]. 

Bagul, Borowy and Salameh [6-8] have developed several 
methodologies for optimally sizing a wind/PV system associated 
with a battery bank for a given load. These methods are based on 
the use of long term data for both irradiance and wind speed. 
However, such approaches are penalized by CPU time due to wide 
data range. Several studies have used the average hourly wind 
speed data over a few years simulation period, but this vision 
strongly filters wind turbine powers. Other researchers [13,14] 
have developed probabilistic methods to determine the annual 
energy of a wind system. In particular in [15-19], authors have 
selected the optimal combination and sizing of wind generators, 
PV modules and storage batteries. 

This paper suggests a systemic methodology for designing the 
optimal combination and sizing of passive wind turbines asso- 
ciated with electrochemical storage. Generally, deterministic opti- 
mization approaches neglect the effects of environmental inputs 
uncertainties (including variation or perturbation of wind speed 
and/or load profiles) which can lead to drastic change of optimal 
solution quality. Several studies [20-26] have taken into account 


Neels number of cells in series 

Neu number of the battery bank changes over 20 years 
Neyc equivalent number of cycles 

Nyc equivalent number of cycles 

Pig extracted power (W) 

Pioad load power (W) 

Pwr wind turbine nominal power (W) 
SOC battery cell state of charge 

Tai battery bank time autonomy (h) 
TCO total cost of ownership (k€) 

Vo battery element nominal voltage (V) 
Vi23 generated wind speed (m/s) 

Vref reference compact wind speed (m/s) 
Virina wind speed (m/s) 

Weyc weight of a cycle 

T repeated wind cycle (days) 

Top operating period (days) 


the increased need for sensitivity analyses to perform a robust 
system design in various research areas.Thus, this study is parti- 
cularly focused on the sensitivity analysis of a set of “good” 
solutions to obtain an optimum system quite insensitive to 
environmental variations. 

The remaining of the paper is organized as follows. The passive 
wind system and the battery characteristics are described in 
Section 2. In the third section, the systemic optimization proce- 
dure is presented. Section 4 is dedicated to the local optimization 
procedure. The models and the systemic optimization formulation 
are exposed in Section 5. Section 6 is reserved to the results and 
the sensitivity analysis. 


2. Description of the system 


The considered system is a “low cost” full passive wind turbine 
(WT) battery charger (Fig. 1) without active control and with 
minimum number of sensors as studied in [27,28]: this local 
optimization loop is not detailed in this paper but only referred 
to previous studies. An experimental prototype of the optimized 
system, especially the PM synchronous generator, has been built in 
LAPLACE Lab and had confirmed the ability of the passive structure 
to draw a power close to the optimal range with small losses. This 
prototype and subsequent study is detailed in [29]. The wind 
turbine parameters have been obtained by applying similitude 
relationships with reference to a 1.7 kW wind turbine which had 
been previously optimized by the local IOD process in [30,31]. 

Both wind speed and load profiles used in this study are 
considered as deterministic data. These data were acquired at a 
typical farm in Tunisia. The load profile is set on 24h and day by 
day repeated (Fig. 2). The wind speed profile has been obtained 
from a previous study [32] by applying a “compact synthesis 
process” on an actual wind speed profile of 200 days duration 
considered as reference data, with the aim of generating a 
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Fig. 1. WT system with battery for stand alone application (rural site 
electrification). 
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Fig. 2. Typical farm load profile for one day. 


25 


20 


Vind (m/s) 


0 1 L L L 
0 2 4 6 8 10 


Time (days) 


Fig. 3. Wind speed compact profile. 


Table 1 
Basic characteristics of a Yuasa NP 38-12I lead acid battery element. 


Nominal capacity C3 30.3 (Ah) 
Nominal voltage Vo 12 (V) 
Nominal discharge current [3 10.1 (A) 


compact profile on a reduced duration of 10 days (Fig. 3) to 
accelerate the optimization process. It should be noted that this 
process synthesizes a wind cycle which fulfils the main character- 
istics of the reference cycle on 200 days, especially the maximum 
wind speed, the average cubic wind speed and the energy content 
of the wind cycle [32]. 

In this study, a lead acid Yuasa NP 38-12I [33] is considered as 
battery element. The basic characteristics are summarized in 
Table 1. 

NB: in the following, let us note that the 12 V Yuasa “battery 
element” is constituted of 6 “cells” of 2 V in series. The association 
of several “battery elements” (12 V) in series and in parallel will be 
named as “battery bank”. 


3. Systemic optimization procedure 


In order to perform the optimization of the whole system 
including the passive WT and the storage battery considering wind 
speed and consumption profiles, we have adopted an approach 
based on two optimization levels: 


Level 1 Systemic Optimization (SO) of the whole system. 
Level 2 Local Optimization (LO) of the WT device. 


Given a wind speed profile, the SO approach consists in 
selecting a WT in the range normally provided by the WT 
manufacturer (WT T,-T, and the corresponding Permanent Mag- 
net Synchronous Generator (PMSG) G,-G,: in our case, these 
devices are obtained from the LO process as presented in [31]). 
The couple (T; G;) and the corresponding storage sizing have to 
satisfy a given load demand at lowest TCO. This compromise can 


i Optimization i : 
: (LO) it 


Systemic 
: Optimization 
(SO) 


Fig. 4. Systemic Optimization process. 


WT and PMSG 
Sizing by similitude 


Fig. 5. Passive WT synoptic. 


be obtained by solving the optimization problem illustrated in 
Fig. 4. 


4. Local optimization 


The aim of this second optimization level is to build a range of n 
generators corresponding to n WTs which will be used in the first 
optimization level (systemic approach). In this approach, the 
battery bank voltage is as constant and equal to 48 V (i.e. a series 
connection of 4 Yuasa 12 V elementary batteries) (Fig. 5). 

The “mixed reduced model” considered in this optimization 
process is detailed in [27,31]. 

In [31], an IOD method, based on a multiobjective optimization, 
has been developed for sizing the elements of a 1.7 kW passive 
wind turbine system. The WT range with related generator 
parameters for various nominal powers has been obtained by 
applying similitude relationships with reference to the 1.7 kW 
wind turbine system [30]. Fig. 6 shows the extracted powers Pex: of 
the passive wind turbines (till 16 kW) obtained by similitude from 
the optimized reference passive wind turbine of 1.7 kW 
(Appendix). It can be seen that the quality of wind power 
extraction of these passive configurations (solid curves) matches 
very closely the behavior of active wind turbine systems operating 
at optimal wind powers by using an MPPT control device (i.e. the 
cubic curves in Fig. 6). It should be noted that the selected wind 
turbine range (1.7 up to 16 kW) is constrained by the average load 
power (Piggq), Since we have adopted a range between % (Piggq) and 
5 (Pioaa). Furthermore, the use of similitude models from a 
reference wind turbine device of 1.7 kW limits the range of power 
expansion to fulfill validity domain of similitude relationships. 
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Fig. 6. Extracted power of the range of passive WT systems obtained by similitude from the reference 1.7 kW WT. 


5. Systemic optimization 


The aim of the SO stage is the minimization of the TCO on a life 
cycle of 20 years of the passive WT coupled with a storage bank 
ensuring the electrification of the isolated farm under a specific 
compact wind speed cycle [32]. To achieve this optimization 
process, a cost model of each system component has to be 
determined. 


5.1. WT cost model 

Generally, the WT subsystem is composed of a turbine, a 
nacelle, a tower, a rectifier, a mechanical transmission system. 
There is no single component that dominates the WT cost. Typical 
owning costs given by “eaglewestwind” [34] for a range of turbines 
from 2 kW up to 20 kW are shown in Fig. 7: system costs include 
costs due to all components of the device. A linear interpolation of 
these costs has led to the following WT cost model: 


Cwr = 1.7Pwr +3 (1) 


The lifetime of WTs being assumed to be at least 20 year, the 
owning cost of WTs will simply be due to the investment cost, 
which assumes that no maintenance costs are needed due to the 
robustness of the full passive wind turbine structure. Such 
assumption can be justified for low power wind turbines, espe- 
cially as the passive structure is particularly robust. 


40 


Cwr=17Pwr+3 


Cwr (k€) 


20 
Pwr KW) 


Fig. 7. Assumption of a WT cost model versus rated power Pwr. 


5.2. Battery bank cost model based on cycling effects 


The battery cost depends on the charge and discharge cycling 
in terms of number and cycle depth that the battery can provide 
(battery lifetime strongly depend on both depth and rate of 
discharge). The equivalent number of cycles (N.yc) is obtained by 
a weighted sum of the cycles in terms of their depth of discharge 
(DOD). The weights are derived from the manufacturer curve 
giving the number of complete cycles that the storage element is 
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capable of delivering in terms of its DOD (characteristic called 
“cycle to failure”). 

The lifetime model, for Yuasa NP 38-12I lead acid battery, uses 
an exponential curve fitting based on the available cycle to failure 
(Cr) curve (Fig. 8) versus DOD [35,36,37]: 


Cr = 177.77 +7807.39.e7873-P0D 2 


Considering the number of cycles provided by the battery with 
a DOD of 100% as a reference, the weight of a cycle depth DOD 
(Wcyc(DOD)) is expressed by the following equation: 


_ C,(100%) 
Weyc(DOD) = ‘C(DOD) (3) 
Thus, the equivalent number of cycles is obtained by: 
Neye = $ Weyc(DOD) x N(DOD) (4) 
DOD 


where N(DOD) is the number of cycles at a given DOD obtained by 
the Rainflow method. 

For an operating period top of 20 years and a repeated wind 

cycle z of 10 days, the approximate cost of the battery bank, taking 
account of cycling constraints over 20 years, is expressed by the 
following equation: 
Cpar(k€) = Neet_sNcet_pCo x 10? Nie (5) 
where Nce;_s is the number of battery cells of 2 V associated in 
series (here considered as constant and equal to 24 to constitute a 
battery of 48 V), Nce_p is the number of battery cells associated in 
parallel and Co is the cost of one deep battery cycle (DOD= 100%) 
estimated at 0.1 € [33]. This latter cost comes from the Yuasa 12 V 
battery element cost which is 108 €. Then, the 2 V battery cell cost 
is Cceu= 108/6 € for 180 allocated deep cycles. Finally, one cell deep 
cycle cost is Ccey/180=0.1 €. 

It should also be noted that the initial purchase cost of the 
battery Cgaro is defined as following: 


Cparo(k€) = Neet_sNcet_pCcelt x 10°? (6) 


5.3. Optimization problem formulation 


The problem is to develop a systemic approach that allows 
optimizing system configurations (passive WT with storage bank) 
that satisfy load power demand with minimum TCO. The TCO 
calculation includes the WT owning cost and battery bank owning 
cost. The systemic optimization process is detailed in Fig. 9. By 
means of the compacted wind cycle obtained from the synthesis 
process proposed in [32] (as presented in Fig. 3), only 10 days of 
system operation have to be simulated in this optimization 
process. Results obtained from this typical but representative 
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Fig. 8. Cycle to failure curve for Yuasa accumulator cell. 
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Fig. 9. Systemic process of optimal design. 


compact profile are then extrapolated over 20 years in order to 
estimate the TCO due to cycling effects. 


5.3.1. TCO 

A single objective function is used consisting in the TCO which 
includes WT and battery bank owning costs over a duration period 
of 20 years. 


TCO = Cwr + Cear (7) 


5.3.2. Design variables 
The optimization problem only uses two discrete design variables: 


© Ncel_p: The number of battery cells associated in parallel (Nees 
being the number of battery cells (2 V) in series is fixed at 24 
for a 48 V battery bank). 

© index i: This index identifies the couples of turbines and associated 
generators (T; and G;) used in the simulation model for computing 
the system TCO and the optimization constraints. Note also that 
this index corresponds with power level (in kW) of the WT. 


5.3.3. Optimization constraints 


è g: is a constraint related to the maximum discharge current 


(Idisch_max)- 
The battery cell maximum discharge current (max(Ice;)), Must 
be less than the current Igisch_max: 


&= max(I¢e1)—Ldisch_ max S 0 (8) 


© g: is a constraint related to maximum charge current (Ich_max)- 
The absolute value of the maximum battery cell charge current 
( min(Icer) | ) must be less than the current Icn_max! 


§2= min(Icer) len max S 0 (9) 


Note that charge/discharge current are respectively considered 
as negative/positive. 
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è g3: is a constraint related to the battery cell state of charge 


(SOC). 


Note that Ncy is the number of the battery bank changes 
necessary for 20 years. Ncp can be calculated from (5) and (6) by 
deriving the following equation: 


The minimum value of the battery cell state of charge min(SOC 


(t)) must be greater than 20%: 
£3 =20%— min(SOC(t)) < 0 


In this study, Igisch_max and I-h_max have been chosen equal to the 


nominal discharge current J3=10.1 A (Table 1). 


As defined, the optimization problem is sufficiently simplified 
(one objective, two design variables and three constraints) to 
justify the use of an exhaustive search (80.14= 1120 combinations) 
instead of adopting more sophisticated optimization algorithms. 
Other approaches may be used in more complex case study, as for 
hybrid systems with several sources. Note also that this simplified 
and efficient optimization problem is mainly due to the structure 
of the design process with 2 loops, one for the WT device and 


another one for systemic approach. 


6. Results and sensitivity analysis 


6.1. Optimization results 


Fig. 10 shows a set of solutions resulting from the systemic 
optimization process, corresponding to a TCO lower than 115.5 k€. 
It should be observed that the TCO range is displayed with a 
reduced scale (1/115 k€) so that TCOs of these solutions can be 
considered as quasi-equivalent. Each point corresponds to one 
solution, (i.e. one couple of variables: index i representing the WT 
rated power and Neer-p). The cheapest solution (circled in Fig. 10) 
corresponds to a configuration with 63 branches in parallel (i.e. 
252 “12 V Yuasa battery elements” for the battery bank) with a WT 
of 13 kW of nominal power (i=13). Table 2 gives the different 
characteristics of some quasi-equivalent optimal solutions and 
Fig. 11 shows the state of charge SOC and the battery cell current 
(cet) evolutions of the solution Sol 13əpt. One can see that a 
reduced DOD range is obtained from that optimal sizing taking 


account of cycling effect on system costs. 


Table 2 
Quasi-equivalent optimal solutions characteristics. 

Sol Sol Sol Sol Sol Sol 
16(a) 16(b) 12 13(a) 13opr_~—- 15a) 

Neel p 49 59 61 62 63 69 
Cwr (k€) 30 30 23 25 25 28.5 

Ncyc (on 20 years) 714 594 622 602 589 520 
Cear (K€) 84.1 84.2 91 89.7 89.1 86.2 
TCO (k€) 114.3 1144 1145 1148 114.2 114.7 

Cpato (K€) 21 25.4 26.3 26.7 27.2 29 

Ncu (on 20 years) 3 3 3 3 3 2 

115 r r 
Sol 13(a) 
Sol 15(a) J 
114.8 k 
z 1146+ ea J 
S gal soie = PO? Sol 15(b) 
Puls 
114.2 } J 
114 , fl , , , , 
45 50 55 60 65 70 75 80 
Neel-p 


Fig. 10. Similar solutions with equivalent TCO. 
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S gar Over 20 rean] 
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(10) Csaro 

Another criterion necessary for the optimal solutions analysis is 
the battery bank autonomy time Taut of each solution. This time is 
defined as the corresponding duration of a complete battery 
discharge up to 20% of SOC, considering the load profile without 
any wind power production. Table 3 displays the autonomy time of 
each solution. 

Fig. 12 confirms the example of Fig. 11 by illustrating a major 
trend: optimization of TCO strongly limits variations of battery 
SOC (reduced values of DOD) in order to enhance its lifetime. 
Results would have been completely different if the cycling was 
not taken into account through the system TCO. In such a case, the 
whole range of DOD would be exploited with deep discharge 
cycles. 

The analysis of the results summarized in Table 2 illustrates 
that some optimized solutions with different characteristics have 
almost the same TCO. The solutions indexed 15(a), 15(b) 
(Pwr=15 kW) and 16(a), 16(b) (Pwr= 16 kW) correspond with the 
cheapest battery owning cost Cgar among others configurations. 
The choice of the solution indexed 13 is optimal in terms of TCO, 
but nearly equivalent and surely not the most relevant solution 
regarding other criteria as the system autonomy. Fig. 12 shows the 
curves of the probability density functions (pdf) of the number of 
cycles versus the depth of discharge (DOD) and their locations on 
the corresponding “cycle to failure” curve for each solution. The 
circle diameter on the “cycles to failure” curve is proportional to 
the cycle occurrence, i.e. the equivalent pdf. These figures show 
that the maximum of DOD for all solutions does not exceed 23% of 
the total storage capacity of the battery bank. These results explain 
why there is a trend to increase the number of battery cells in the 
solution Sol 15(a) compared with the solution Sol 16(a): this is 
obviously to further expand the total storage capacity of the 
battery bank and to reduce subsequently deep discharge cycles 
therefore enhancing the lifetime. 
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Fig. 11. SOC and Ie; evolutions for the optimal solution. 


Table 3 
Time autonomy of quasi-equivalent optimal solutions. 


Sol 16(a) Sol 16(b) Sol 12 Sol 13(a) Sol 130pt Sol 15(a) Sol 15(b) 


Tau (h) 5 6 62 64 6.6 71 8.1 
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Fig. 12. pdf of Ncyc versus DOD on cycle to failure curve of each solution. 
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Table 4 
Wind speed profiles characteristics. 
c k Vmin Vmax (V) Vv’) 
Vref 95 2.3 01 25 8.3 1195 
Vi 9 2 0.04 25.3 8 867 
V2 10 2 0.04 31.6 8.9 1189 
V3 12 2 0.1 34.5 10.8 2072 


Furthermore, at equivalent system cost, it is certainly prefer- 
able to dispose of a greater storage capacity in order to ensure load 
demands during longer periods without wind (i.e. autonomy 
Table 3). More generally, it is interesting to analyze these solutions 
by sensitivity analysis versus changes of environmental data (wind 
speed and load) then to rebuild the optimization problem to 
choose a robust and optimal solution. 


6.2. Sensitivity analysis of optimal solutions versus environmental 
data variations 


This section is devoted to the sensitivity analysis of the 
equivalent optimal solutions presented in Fig. 10 versus environ- 
mental input changes (i.e. wind speed and load profile variations) 
with the aim of identifying the most robust solutions. 


6.2.1. Sensitivity versus wind speed profile changes 
Equivalent optimal solutions represented in Fig. 10 have been 
obtained for a particular wind speed compact profile (V,ef) (see 
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Fig. 13. Probability distribution functions of the wind speed profiles. 


Fig. 3). Vref distribution corresponds to the characteristics of a 
Weibull law with a scale factor c=9.5 m/s and a shape factor 
k=2.3. Here, we aim at analyzing solution robustness versus 
statistic characteristics mentioned in Table 4. 

For this purpose, three particular time profiles of wind speed with 
200 days duration are generated from different Weibull distributions 
(Vi, V2, V3) as displayed in Fig. 13 and characterized in Table 4. 

Table 5 gives the number of cycles over 20 years Ncyc, the 
battery bank cost Cgar and the TCO of several equivalent optimal 
solutions for each wind speed profile. These results show that the 
solution Sol 15(b) is the most robust in terms of TCO change: (Sol 
15(b) remains lower than the one obtained with other solutions 
whatever the wind conditions. This robustness is partly due to the 
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Table 5 


Wind speed profile variations results. 


Sol 12 


Vref Neye 714 594 622 


114.38 


TCO (k€ 107 94 104 


Sol 16(a) Sol 16(b) 
Cear (KE 84.1 84.2 91.1 
TCO (k€ 114.32 114.40 
Vı Neye 1015 829 995 
Cear (KE 119 117 145 
TCO (k€ 149 147 169 
V2 Neye 846 659 806 
Cpar (K€ 99 93 118 
TCO (k€ 129 123 141 
V3 Noyce 653 455 557 
Cpar (K€ 76 64 81 
Table 6 
Wind speed profile V characteristics. 
c k Vmin Vmax (V) (Vy 
V 9.5 2.3 0.18 25.67 8.41 914 
8 T 
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Fig. 14. Load profile changes. 


battery capacity of the Sol 15 (b) (Nee p=79) which is the highest 
among all other solutions. 


6.2.2. Sensitivity versus load profile changes 

The sensitivity analysis of the optimal solutions versus load 
profile is performed by considering an actual wind profile of one 
year duration (V) whose characteristics are given in Table 6. 
Changes made with respect to the original load profile consist in 
inserting a variation of 25% of load power during the even days 
and —25% during odd days (Fig. 14). 

Table 7 shows that the Sol 15(b) remains the most robust 
solution in terms of TCO versus variations of load profile. As 
mentioned in the case of wind speed profile variations, this is 
explained by the highest storage capacity. 

Beyond this sensitivity analysis, we can guess that the solution 
with the highest battery capacity (Sol 15 (b)) is also the most 
autonomous in case of wind drops: it is subsequently the most 
reliable and relevant for this application. 


7. Conclusion 


This paper illustrates a systemic optimization approach 
devoted to the optimal design of a full passive WT with storage 


Sol 13(a) Sol 13 opt Sol 15(a) Sol 15(b) 
602.98 589 520 452 
89.7 89 86.2 85.8 
114.86 114.25 114.74 114.37 
924 900 728 623 
137 136 120 118 
162 161 149 146 
752 736 582 480 
111 111 96 91 
137 136 124 119 
499 492 402 340 
74 74 66 64 
99 99 95 92 
Table 7 
Sensitivity versus load profile variations. 
Sol 16 Sol16 Sol Sol 13 Sol Sol 15 Sol 
(a) (b) 12 (a) 13op (a) 15(b) 
Original Neyc 824 655 842 761 746 578 495 
profile Cear 97 92 123 113 112 95 94 


(k€) 
TCO 127 122 146 
(k€) 


138,4 138 124 122 


Modified Neye 791 650 811 746 728 573 479 

profile Caii 93 92 118 111 110 94 90 
(k€) 

TCO 123 122 142 136 135 123 119 


(k€) 


dedicated to stand alone applications (here for rural electrification 
purpose). System environment, especially wind and load condi- 
tions is integrated in the systemic design by optimization. The 
integrated optimal design problem was divided into two optimiza- 
tion processes: a local optimization which aims at designing a set 
of optimized WT in order to constitute a manufacturer range to be 
used in the second optimization process dedicated to the whole 
system design. The systemic optimization objective is to minimize 
the total owning system cost, integrating simultaneously WT and 
battery costs by taking account of the number of device change 
due to cycling effects. It is interesting to note that such a systemic 
vision of the TCO over the whole life cycle (here 20 years) leads to 
oversize batteries in order to reduce DODs (lower than 20%) in 
order to reduce the number of device changes over the 20 years. 
Optimal results were analyzed and discussed then completed by a 
sensitivity analysis which has compared several selected solutions 
previously considered as quasi-equivalent. The sensitivity of these 
quasi-optimal solutions has been analyzed versus changes of 
environmental data (wind speed and load profiles variations) 
showing that solutions with highest storage capacity were the 
most relevant, being the most robust and also the most autono- 
mous in case of lack of wind. The availability of the developed 
methodology was successfully demonstrated with the considered 
environmental data and it is appropriated for other acquired sets 
of experiment data to find consistently the optimal generation 
system design. 
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Table A1 
Turbine: rated power Radius (m) Friction coefficient (N m s/rd) 
Reference: 1.7 kW 1.25 0.025 
Ti: 1 0.95 0.0087 
Ta: 2 1:35 0.034 
Ts: 3 1.66 0.07 
Ta: 4 1.91 0.13 
Ts: 5 2.14 0.21 
Ts: 6 2.34 0.31 
Tz: 7 2.53 0.42 
Tg: 8 2.71 0.55 
Tg: 9 2.87 0.70 
Tio: 10 3.03 0.86 
Ty: 11 317 1.04 
Tiz: 12 332 1.24 
Tiz: 13 3.45 1.46 
Tia: 14 3.58 1.69 
Ty5: 15 3.71 1.94 
Tig: 16 3.83 2.21 
Table A2 
PMSG: rated power (kW) Flux (Wb) Resistance (mQ) Inductance (mH) 
Reference: 1.7 0.20 12626 1.34 
Gi: 1 0.15 2784 1.57 
Gz: 2 0.21 94.3 1.36 
G3: 3 0.28 60.9 153 
G4: 4 0.33 40.4 15 
Gs: 5 0.35 26.5 1.33 
Ge: 6 0.42 24.2 1.55 
G7: 7 0.41 15.6 122 
Gg: 8 0.47 14.6 1.37 
Go: 9 0.53 13:7 1:51 
Gio: 10 0.47 8.3 1.05 
Gun: 11 0.51 7.9 1.14. 
Gia: 12 0.56 7.6 122. 
Giz: 13 0.61 73 13 
Gia: 14 0.65 7 1.38 
Gis: 15 0.7 6.8 1.46 
Gis: 16 0.75 6.6 1.54 
Appendix 


See Tables A1 and A2. 
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